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Abstract 

The probablity current for a quantum spinless relativistic particle is introduced based on the 
Hamiltonian dynamics approach utilizing the Salpeter equation as an alternative of the Klein- 
Gordon equation. The correctness of the presented formalism is illustrated by examples of exact 
solutions to the Salpeter equation including the new ones introduced in this work. 
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I. INTRODUCTION 



The problem which arised in the very early days of quantum mechanics was finding 
the relativistic counterpart of the Schrodinger equation. The most popular choices of the 
relativistic quantum mechanical equations for spinless massive particle are the Klein-Gordon 
equation [IHl] and less familiar relativistic Schrodinger equation which is usually referred to 
as the spinless Salpeter equation [5ti23] . The latter can be regarded as the "square root" of 
the Klein-Gordon equation and is based on the approach which is sometimes referred to as 
the relativistic Hamiltonian dynamics [211 [25] . 

The advantage of the Klein-Gordon equation is that it is manifestly covariant. Its well- 
known fiaw related to the fact that this equation is of second order in the time derivative, 
is the problem with the probabilistic interpretation. Namely the time component of the 
probability four-current which in the nonrelativistic case should coincide with the probability 
density, can be negative. As a matter of fact, as the second order equation the Klein-Gordon 
equation can be cast into system of first order equations which can be represented in matrix- 
form equation. An example is the Kemmer equation [261 EZj and the two-dimensional system 
discussed by Feshbach and Villars ^28j. Nevertheless, the reduction of this kind does not 
solve problems with the Klein-Gordon equation and generate new ones. Furthermore, the 
appearance of both signs of energies for the solutions to the Klein-Gordon equation leads 
to the occurence of the "Zitterbewegung" and the Klein paradox. The " Zitterbewegung" 
or trembling motion i.e. rapidly oscillatory motion whose amplitude and period are of order 
h/mc and h/mc^, respectively, is the result of the interference between positive and negative 
energy states. The Klein paradox relies on possibility of transmission to states with negative 
kinetic energy in the electrostatic step-like potentials. To circumvent the problem one usually 
suggests the creation of the particle-antiparticle pairs, that is the need is indicated for second 
quantization as with the Pauli-Weisskopf approach [25] where the probability current is 
reinterpreted as a charge current. Finally, the scalar product in the space of solutions of 
the Klein-Gordon equation, although Lorentz covariant, is not positive definite which causes 
serious problems. Therefore, the severe difficulties arise with the physical interpretation of 
the Klein-Gordon equation and it is widely believed that based on this equation one cannot 
construct consistent one particle relativistic quantum mechanics. 

The disadvantage of the Salpeter equation is that it is not manifestly covariant. Another 
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problem which is sometimes indicated is the nonlocahty of the relativistic Hamiltonian 
which is the pseudodifferential operator. However, as was pointed out by Lammerzahl [22] 
the nonlocahty of the Salpeter equation does not disturb the light cone structure. The 
macrocausality of the second quantized version of this theory was also reported therein. In 
addition, it has been demonstrated by Foldy [6J that the space L'^{M.^ ,d^x) of solutions to 
the Salpeter equation is invariant under the Lorentz group transformations which is a direct 
consequence of a unitary relationship between L'^{M.^,d^x) and the space of the unitary 
irreducible representation of the inhomogeneous Lorentz group [25J. On the other hand, 
it is clear that nonlocal pseudodifferential equations like the Salpeter equation are much 
more complicated than the local differential ones like the Klein-Gordon equation. In spite 
of difficulties, the great advantage of the Salpeter equation is that it possesses solutions of 
positive energies only, so we have no problems with paradoxes mentioned above occurring in 
the case of the Klein-Gordon equation. We also point out that agreement of predictions of 
the spinless Salpeter equation with the experimental spectrum of mesonic atoms is as good 
as in the case of the Klein-Gordon equation [30] • Moreover, the possibility of probablistic 
interpretation in the quantum case as well as clear classical physical content of the Salpeter 
equation was the motivation for its wide usage in the phenomenological description of the 
quark-antiquark-gluon system as a hadron model [3T| [32]. 

In this work, after discussion of the relationship of the Salpeter equation with the cor- 
responding integro-differential equation, we introduce and examine the probablility current 
derived from the spinless Salpeter equation. We show that such current has all good proper- 
ties of its nonrelativistic counterpart. In particular, our analysis shows that the nonlocahty 
of the Salpeter equation does not disturb causal propagation of particles. The theory is 
illustrated by the concrete examples of exact solutions of the Salpeter equation. 



II. THE SALPETER EQUATION 

This section is devoted to the discussion of the basic facts about the Salpeter equation. 
The Hamiltonian of a relativistic classical particle subject to the potential V{x) is 

H = y^c'^p^ + m'^c^ + V{x), (2.1) 
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where m is the mass of the paricle, c is the speed of hght and x is the three-position. In 
relativistic quantum mechanics the system defined by the Hamiltonian (2.1) is described by 
the Salpeter equation of the form 



.^^^v^g^ = [Vm^c^ - n^c^A + V{x)]cj){x, t), (2.2) 

where A = V^. The Eq. (2.2) is obtained by the quantization procedure utihzing the 
Newton- Wigner locahzation scheme [6]. In this scheme we have the standard quantization 
rule X ^ X, and p — )■ — i/iV. As a consequence, the space of solutions of the Salpeter 
equation is the Hilbert space L'^{M.^ ,d^x) with the scalar product 

{^\iP) = j £x(P*{x)i;{x). (2.3) 

Therefore, according to the quantum-mechanical spirit, we should identify \(f){x,t)\'^ with 
the probability density p{x,t) satisfying the normalization condition: 

d^xp{x,t) = l. (2.4) 

Motivated by usage of limiting procedures in the nonrelativistic and ultrarelativistic case 
as well as some dimensional considerations we keep in this section and the following one 
the physical constants h and c. The natural units h = 1 and c = 1 are utilized in Sec. IV. 
Performing the Fourier transformation 

cf>{x,t) = — ^3- / d'pe''^4>{p,t), (2.5) 
(27r) 2 h'^ J 

(in the following, whenever it is clear from the context, we omit designation of a region of 
integration) we obtain from (2.2) the following equation 

.^90^ = Wm^c^+P^c^ + ViihVpMp, t). (2.6) 

It is clear that (2.6) is the partial differential equation of finite order only for V{x) polynomial 
in X. Note that a/^tj^cM-~P^ is the so called "symbol" of the psudodifferential operator 
y/mP'c'^ — h^c^A, that is 

Vm^c^ - h^c^A(P{x,t) = ^ [ d^p^rn^c^ +p'^c^e''^4){p,t). (2.7) 

(27r) 2 h? J 

On taking the inverse Fourier transformation 



and making use of the identity [33] 

Vx2 + a? e'^'" dx = ---KMp\), (2-9) 

m 

where Kij{z) is the modified Bessel function (Macdonald function), as well as the differenti- 
ation formula satisfied by the Bessel functions such that 

K[{z) = -^K^{z)-K2{z), (2.10) 

we get from (2.7) the following formula on the action of the pseudodifferential operator 
y/rrL^d^ — h^c^A in the coordinate representation: 

V^^^,i.,.).l,'yKi.-yM (2.11) 

where the function K{x — y) is 

and \a\ designates the norm of the vector a. Thus, it turns out that the pseudodifferential 
operator Vm^c^ — h?c^A can be defined as the integral operator with the kernel (2.12). 
Consequently, the Salpeter equation (2.2) takes the form of the integro-differential equation 

.^d^i^ = I d'y K{x - y)<P{y, t) + V{x)<P{x, t). (2.13) 

We also remark that the nonlocality of the Salpeter equation is related only to the kinetic 
energy term described by the integral operator from the right-hand side of Eq. (2.13) and 
does not depend on the potential. Therefore, the nonlocality is not connected with potential 
forces acting on a quantum particle. Consider now the massless limit m = 0, when the 
Salpeter equation is 

.d(f){x, t) 



ih—^-^ = [hcV^ + V{x)](j){x, t). (2.14) 
Of course, the corresponding Fourier transform (l){p,t) fulfils the massless limit of Eq. (2.6), 



I.e. 



= + ViihV^mp, t). (2.15) 



dt 

Taking the limit m — )■ of Eq. (2.12) and using the asymptotic formula 



K2iz) = ^, z^O, (2.16) 



we find for m = 



K{x-y) = ^ 



14' 



(m = 0). 



TT^ \x — y\ 

Finally, consider the simplest case of a relativistic massless particle moving in 
the Salpeter equation is 



d4>{x, t) 
dt 



(f){x,t). 



On performing the Fourier transform 



(j){x,t) = 



we obtain the counterpart of Eq. (2.6) 



27rh J- 



ipx 



at 



d_ 
dp 



Clearly, the one-dimensional version of (2.7) is 

~W . 1 



m'^c'^ — h?c'^—^ (bixA) = 



dp a/ im?d^ + p^c^ e ^ 0(]9, t). 



Taking the inverse Fourier transformation 



(l>(P^t)^^= / dye (l){y,t), 



and using (2.9) we find 

\jm?d^ - fPc^—(f){x,t) = J dyK{x - y)(f){y,t), 

where the kernel K[x — y) is given by 

K[x -y) = -. tAi [^\x-y\] . 

n \x — y\ \ n / 



y\ 

It follows immediately from (2.24) and the asymptotic formula 

Ki{z) = -, z^O 
z 

that the kernel in the massless case is of the form 

ch 1 



K{x - y) 



The relation (2.26) is also a direct consequence of (2.21) with m = 0, (2.22), and the identity 



As far as we are aware, the first example of the exact solution to the Salpeter equation, 
referring to the massless free particle moving in a line, was considered by Rosenstein and 
Horwitz |TT] whereas the solution for the massive free particle was discussed by Rosenstein 
and Usher [12] . The WKB approximation technique to the Salpeter equation was developed 
in Refs. [8] and f9]. In a series of papers by Hall, Lucha and Schoberl (see for instance Ref. 
[T5]). the energy bounds for the Salpeter equation were analyzed. To our best knowledge, the 
first example of the nontrivial exact solution to the Salpeter equation in the case of a particle 
in M^, was our solution to the Salpeter equation for a relativistic massless harmonic oscillator 
[23] . More precisely, we derived the exact stationary wavefunctions and corresponding exact 
spectrum of the energy expressed by means of zeros of the Airy function. The correctness 
of the quantization based on the massless Salpeter equation was confirmed by the good 
behavior of the related probability density and expectation values of observables. Recently, 
the new exact solution to the Salpeter equation has been reported referring to the free 
massive particle on a line, with the Gaussian initial wavefunction [S^ . This solution is very 
complicated and it has the form of infinite power series expansion with coefficients expressed 
by means of integrals of special functions. In Sec. IV we introduce other new examples of 
exact solutions to the Salpeter equation. In particular, we derive the nontrivial solution of 
this equation for a massless particle in a linear potential. 

III. PROBABILITY CURRENT 

In this section we introduce the probability current for a quantum spinless relativistic 
particle and discuss its basic properties. We first discuss the probability density and the 
probability current for the Klein-Gordon equation such that 



where the potential V{x) is introduced in the equation by means of the vector minimal 
coupling scheme. We point out that the denomination "square root" of the Klein-Gordon 
equation, mentioned in the introduction is appropriate for the Salpeter equation only in the 



m 




(2.27) 




(3.1) 
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case of a free particle. Indeed, by squaring (2.2) we obtain 

b{x, t) = ^m^c^ - h^c'A + [Vm^c^ - h^d'A, V{x)] | ^{x, t). (3.2) 



We recall that in the case of the Klein-Gordon equation the probability density is given by 
(compare Ref. [35] ) 

ih f i^dip , dip* ^ 2i^^| , 2 



Since the Klein-Gordon equation is second order in time derivative, therefore the initial 
values of ijj and ^ can be arbitrary. We conclude that pkg can be either positive or negative 
and the problem arises with the probablistic interpretation of the Klein-Gordon equation. 
The expression for pkg reduces to the nonrelativistic form in the nonrelativistic limit c 
oo. However, to show this one should first assume the validity of the Salpeter equation 
(2.2) or equivalently restrict to positive energy solutions of the Klein-Gordon equation only. 
The formula on the probability current for a Klein-Gordon particle is identical with the 
nonrelativistic one, that is we have 

jKG = -;^(^*V^-^V^*). (3.4) 

We stress that neither the probability density (3.3) nor the probability current (3.4) has the 
correct limit m — )■ 0. Furthermore, it can be checked that the continuity equation implied by 
the massless Klein-Gordon equation obtained from (3.1) by setting m = requires existence 
of some universal constant with the dimension of length. This is yet another disadvantage 
of the Klein-Gordon equation. 

We now return to the Salpeter Eq. (2.2). Proceeding analogously as in the case of the 
nonrelativistic Schrodinger equation we find 

+ i(0Vm2c4 - K^c^A (j) - (lyVni^c* - h^c^A (j)*) = 0. (3.5) 

ot h 

Now, the probability density p{x,t) expressed in terms of the Fourier transform (piPjt) can 
be written in the form 

p{x,t) = \^{x,t)\' = J d'pd'ke''-^^*ip,t)4>ik,t). (3.6) 

On using (3.5), (3.6) and the continuity equation 

|h-V.,- = 0. (3.7) 
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where p{x,t) is the probabihty density and j{x,t) is the probabihty current, we arrive at 
the following formula on the probability current: 

It should be noted that the probability current (3.8) has the correct nonrelativistic limit 
c — 7- oo. Namely making use of (2.5) we easily find 

ih 

hm 3 = - — (0*V0 - 0V0*). (3.9) 

c^oo zm 

Furthermore, it follows immediately from (3.8) that 

j{x,t)d'x = {(t)\vcp), (3.10) 



/• 



where v is the operator of the relativistic velocity 



* = (3.11) 



where p = — i/iV, and po = E/c = a/ m^c^ + p^ = \^m?cP~^-h?A. The formula (3.10) is the 
relativistic counterpart of the well-known nonrelativistic expression describing connection of 
the integral of the probability current and average velocity in the given state. We stress 
that the relation (3.10) is not valid in the case of the Klein-Gordon equation. Bearing in 
mind the critique of the Salpeter equation based on its nonlocality, it is also worthwhile to 
point out that the length of the average velocity (3.11) related to the discussed probability 
current via Eq. (3.10), does not exceed the speed of light. Indeed, we have 



|t)0)^ < {<P\v'<P) = c'iM < c^ (3.12) 

p^ + m^c^ 



where the inequality is saturated at m = 0. Yet another good property of the probability 
current (3.8) is the existence of the massless limit. In fact, putting m = in (3.8) we 
immediately get 

j{^,t) = J rfV'fc^^e^^0*(p,t)0(fc,t), (m = 0). (3.13) 

The probability current can be expressed in terms of the solution (f){x,t) to the Salpeter 
equation. Indeed, substituting in (3.8) the Fourier transform (2.8) and using (2.9), (2.10) 
and the identity [33] 



^ a2 - a;2 a 



dx=—sma\p\, (3.14) 



we find 



where 



j{x,t) = j d^y£zK{x;y,z)(P*{y,t)(l){z,t), 



(3.15) 



K{x;y,z) 



(3.16) 



2 S 

im c 



1 



(2'7r)'^/i2 ■ " ■ i\y — x\\x — z\ \y — x\ + \x — z 

Hence, making use of the theorem on the gradient [36j 



d^x'V(p{x) = / dS(p{x), 
V Js 



mc 
LIT 



{\y -x\ + \x- z\) 



(3.17) 



where S is the oriented boundary of the volume V, and taking into account that the Bessel 
functions K^{z) approach zero as |z| — t- oo, we arrive at the relation 

j{x,t) = J d'y£z K{\y - xl\x - z\)[(j)*{y,t)VMz,t) - (f){z,t)Vy(f)*iy,t)], (3.18) 



where 



K{\ul\w\) 



1 

\m c 



mc 
li 



{\u\ + \w\) 



(3.19) 



{2'KYf9 \u\\w\ \u\ + |tf;| 
The formula (3.18) is remarkable. In fact, it means that the relativistic probablility current 
has the form resembling the "smeared" nonrelativistic one. 

We now discuss the case m = 0. On taking the massless limit m — )■ and making use of 
(2.16) we easily get from (3.16) the following formula 

1 1 



K[x;y,z) 



2ic ,^ ^ ^ 



(m = 0). (3.20) 



\y — x\\x — z\ {\y — x\ + \x — z\Y 
The relation (3.20) can be also easily derived from (3.13), (2.8), (3.14), and the identity [33] 

2a 



e{x)sm.ax e*^^ dx 



(3.21) 



where e{x) is the sign function. Proceeding analogously as with (3.15) we get 

2ic 1 1 



K{\ul\w\) 



(m = 0). 



(3.22) 



(27r)3 l**!!'"^! + 

On the other hand, the formula (3.22) is an immediate consequence of (3.19) and the asymp- 
totic formula (2.16). 

Now, whenever the initial wavefunction is real then the solutions to the Salpeter equation 
(2.2) satisfy (p*{x,t) = (f){x, —t). Of course, it is a well-known property of the Schrodinger 
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equation. An immediate consequence of this relation and (3.15) is j{x, —t) = —j{x,t), that 
is j{x,t) is an odd function of time. Furthermore, it follows easily from (3.16) that if the 
wavepacket fulfils (j){—x,t) = 0(x ,t) or (j){—x,t) = — (f){x ,t), then j{—x,t) = —j{x,t). 

We now return to the Salpeter equation (2.2). Using the continuity equation (3.7) and 
the formal power series expansion of the square root from the right-hand side of (2.2) we 
get the following formula on the probability current: 

imc2^(2n-3)!! / h \ 2-1 



J 



h ^ (2n)!! \mc , 

n=l ^ ' ^ ' k=0 



Y.'-^NJr £-J E(-i)'vV'v— V. (3,23) 



We point out that (3.8) can be formally obtained from (3.23) and (2.5). Nevertheless, the 
formula on the probability current (3.23) is mathematically less sound and technically less 
convenient than the integral representation (3.8). Furthermore, in opposition to (3.8), the 
current given by (3.23) has no correct limit for m — 0. 

Finally, let us restrict to the simplest case of a relativistic massless particle on a line. 
Equations (3.6), (3.8) and (3.13) take then the following form 

^ ' dpdke'^-^^*{p,t)4>{k,t), (3.24) 



p{x,t) = 


|0(x,t)|2 


j{x,t) = 




27ih^ /' 


3{x,t) = 







dpdk P + k P^^^d,*(j>A)d,(kA).^ (3.25) 

a/ m?c^ + p2 + V m?c^ + k"^ 

^P^^ ^y^^^^ {^ = ^). (3-26) 

respectively; here we have used the Fourier transformation (2.19). Let us focus our attention 
on Eq. (3.25). Taking into account (2.22), (2.9), and the identity [55] 

00 ipx 

dx = i'Ke{p)e^'''P, (3.27) 

-00 ~^ ^ 

we get from (3.25) 

j{x,t) = J dydzKix;y,z)<P*iy,t)^iz,t), (3.28) 

where 

imc^ e(x — y) — e(x — z) /mc, ,\ , , 

K{x;y,z) = Y -K, (—y-z ). 3.29 

Hence, taking the limit m — > and using (2.25) we find the following formula on the function 
K{x; y, z) for the probability current in the massless case (3.26) 

ic e(x ~ y) — e(x — z) , , , , 

K{x;y,z) = -- _ ^, (m = 0). (3.30) 
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We now discuss the case of the massless particle on a hne in a more detail. Let 0±(a;, t) 
be solution of the Salpeter equation 



for a free massless particle with positive (negative) momentum moving in the line to the 
right (left) . Then we can identify 

0±(p,i) = ^(±p)0(p,i), (3.32) 

where 9{p) is the Heaviside step function. By means of the Fourier transformation (2.19) 
for m = 0, we conclude that 0+ and 0_ depends on. x — ct and x -\- ct, respectively, that is 

0+(x, t) — (f{x — ct), (3.33) 

and 

(l)-{x,t) ^i/jix + ct). (3.34) 

Let furthermore p±{x,t) = \(f>±{x,t)\'^ and j±{x,t) designate the probability density and 
probability current, respectively, corresponding to (f)±{x,t). An immediate consequence of 
(3.24) and (3.26) is 

j± = ±cp±. (3.35) 

We point out that the relation (3.35) is a natural massless relativistic counterpart of the 
well known nonrelativistic formula on the probability curent (see (4.29)) expressed by means 
of the velocity and the probability density. The continuity equation in the one-dimensional 
case such that 

is satisfied identically. Of course, the general solution (f){x, t) to the Salpeter equation (3.31) 
is 

4>{x,t) = -^[M^,t) + M^,t)] = ^b(^ - ct) + iPix + ct)]. (3.37) 

On the other hand, (3.37) is the well-known general solution of the one-dimensional Klein- 
Gordon equation for a massless free particle of the form 



dt'^ dx^ 
12 



0(x,i) = O. (3.38) 



Notice finally, that for massless particle moving to the right or left, the expectation value 
{x{t)) of the position operator coincides with the classical trajectory of a massless particle 
on a line. Indeed, we have 

/oo 
dxxp{x T ct) = ±ct + (x(0)). (3.39) 
-oo 

IV. EXACT SOLUTIONS 

In this section we introduce the exact solutions to the Salpeter equation and discuss 
the corresponding probability density and probability current. We begin with the one- 
dimensional cases. 



A. Free massless particle on a line 

We first study a relativistic free massless particle moving in a line. The corresponding 
Salpeter equation can be written in the form 



i^^^f^.(M), (4.1, 

where we set c = 1. Consider the evolution of the (normalized) wavepacket 

2 at 



This package is referred in Ref. [TT] as to the "Lorentzian" one. On performing the Fourier 
transformation 

0(x,t) = ^ / dpe'P^4>{p,t), (4.3) 



^27r 

where we set h = 1, we get from (4.1) the following equation 

= (4.4) 

subject to the initial condition 

0(p,O) = y^e-'^lfl, (4.5) 
where the use was made of the identity 

oo 1 

e'P^ dx = -e-'^IPl . (4.6) 

X'' + a 
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The solution to the elementary equation (4.4) with the initial condition (4.5) is 

0(p,t) = y^e~ IPl. (4.7) 

Hence, using (4.3) we obtain the normalized wavefunction at any time 

12a a + it , ^ , 

(P{x,t) = \ . — T- 4.8 

' ^ V TT x2 + (a + it)2 ^ ^ 

The solution (4.8) was originally derived in Ref. [H], however the definition of the probability 

current suggested therein is different from ours. An immediate consequence of (4.8) is the 

following formula on the probability density 

The time development of the probability density (4.9) is presented in Fig. 1. Furthermore, 

taking into account (3.26) and (3.36) we find after some calculation 

a (x + t)2 + a2 ax — 3^^ + 

" 4^ (x - ty + a^~^{x^-t^ + a2)2 + 4aH^ " ^ ' ^' 

The time evolution of the probability current (4.10) is shown in Fig. 2. We point out that 

there is no singularity in (4.10) for t = 0. Namely, we have limt^Q j{x,t) = 0. 

We now return to (4.7). In view of (3.32) we have 

^±(p,t) = v^^(±p)e-("+'*)IPl. (4.11) 

Hence, taking into account (4.3) it follows that the wavepacket referring to the particle 

moving to the right and left, respectively, is given by 

Fa ±i , , 

(p±{x,t) = J — . 4.12 

Therefore, the corresponding probability density and probability current are 

p±(x,t) = \Mx,t)\' = ±Mx,t) = ^(^^^j^- (4.13) 

The time development of the probability density p+{x,t) is shown in Fig. 3. 



B. Free massive particle on a line 

We now discuss a relativistic free massive particle moving in a line. The corresponding 
Salpeter equation is 

(4.14, 

14 




-20 -10 10 20 

X 

FIG. 1. The time evolution of the probabihty density (4.9) related to the solution of the Salpeter 
equation for a free massless particle in one dimension. Because of choice of the natural units (c = 1 
and h = 1), the units of p, j, x and t used in the figures referring to the one-dimensional case, are 
m~^, m~^, m and m, respectively. The parameter a = 1 m. 




X 



FIG. 2. The time development of the probability current (4.10) for a free massless particle moving 
in a line. The parameter a = 1 m. 
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FIG. 3. The behavior of the probabihty density given by (4.13), referring to the case of 

the free massles particle moving to the right. The stable maximum of the probabihty density is 
going with the speed of light c=l. 

Performing the Fourier transformation we find 

= ^m^+p^^{p,t). (4.15) 

at 

Consider now the particular solution to (4.15) of the form 



i/2miii(2ma) 



Using the identity [3Z] 

I dx exp(— oa/xM-^) COS7X = ^ - i^i(/3A/a^ + 7^), Rea > 0, Re/3 > 0, 
Jo V + 7^ 

(4.17) 

one can easily derive from (4.16) the following (normalized) solution to the Salpeter equation 
(4.14): 

0(x,t) = I w 2""^'^ K,[m^x^ + {a + m (4.18) 

Y TcKi{2ma) + (a \tf 

Taking into account the asymptotic formula (2.25) we find that the solution (4.18) is a 

generalization of the solution (4.8) obtained for the massless particle to the case m > 0. 

More precisely, we have 

/2a a -|- it , , 

lim X, t) = \ — . 4.19 
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FIG. 4. The time development of the probabiUty density p{x,t) = corresponding to the 

case of the free massive particle, where cj){x,t) is given by (4.18). The mass m = 0.5 and a = 1 m. 

The solution (4.18) was independently derived in Ref. [12j. The probability density p{x,t) = 
|0(x,t)p corresponding to the the wavefunction (4.18) is presented in Fig. 4. It turns out 
that, in opposition to the massless case, the wavefunction 0(x, t) spreads out as time passes. 
Finally, taking into account (3.34), (2.10), the identity 



(4.20) 



and the fact that 0(x,t) is an even function of x (see discussion below Eq. (3.22)) we find 
the following formula on the probability current: 



m 



TiKii2ma) Jq 
a — it 

X 



dx Im 



^x'^ - (a + itf 
x2 + (a + it)2 



K2[Tn^/x^ + {a + ity] - Ko[m^/x'^ + (a + it)^] 



-.Ki[m^/x^ + (a-it)^] } . 



^y x'^ + {a — it) 

The time development of the probablility current (4.21) is shown in Fig. 5. As expected in 
view of the behavior of wavef unctions, the probablility current spreads out. 



(4.21) 



C. Massless particle in a linear potential 

Our purpose now is to analyze the case of a relativistic massless particle moving in a line 
in a linear potential. Classically, this system is defined by the Hamiltonian 



H = c\p\ + nx, 



(4.22) 
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FIG. 5. The plot of the probabihty current (4.21) versus time, where m = 0.5 and a = 1 m. 



where > is a parameter. The corresponding Hamilton equations lead to the solution of 
the form 

<t) = --[\- l^t + p{0)\ -p{0)] + x{0). (4.23) 
A* 

In the following we put c—1. Under simplest initial conditions x{0) — 0, p(0) = we have 
x{t) — —\t\ (a A-shaped trajectory). On the quantum level this system is described by the 
Salpeter equation of the form 



.d(f){x,t) 



Q2 



0(x, t) + ^X(f){x, t), 



dt V dx"^ 

where we set c = 1 and h—1. The Fourier transform of (j){x, t) satisfies 

i^^ = |p|0(p,t) + i/.^^. 
The following solution of Eq. (4.25) can be derived easily 

(j){p,t) = e^^x{P + fJ't), 
where x(p) is an arbitrary function. Now, we choose the initial wavepacket so that 



(4.24) 



(4.25) 



(4.26) 



(4.27) 



where A > is a parameter and C is a normalization constant. Hence, the normalized 
wavefunction in the momentum representation is 



Hp, t) ^ \ — I e' 2m e 2m 



(4.28) 
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Eqs. (4.28) and (4.3) taken together yield the normahzed wavefunction such that 
1 /A/i 



(f){x,t) 



+ 



e 2 



p(At-ij:)^ 



M(-At+i3:)^ 

=e 2{i+i) erfc 



2(1 + i) 



-At + ix) 



(4.29) 



_g 2(l-i) QYfc 



— (At - ix) 



2(1 -i) 

where erfc(2;) = 1 — ^ e~^^dt is the complementary error function. A remarkable prop- 
erty of the corresponding probability density p{x,t) = |0(a;,t)p presented in Fig. 6, is the 
behavior of its maxima following the classical A-shaped trajectory. Moreover, the expecta- 
tion value of the position operator, calculated easily in the momentum representation is of 
the form: 



m) 



TT 



terf ( a/ A/it) 



(4.30) 



where erf(t) = ^ Jq & dr is the error function. Hence, we find the asymptotic formula 



{x{t)) = -|t| 



|t| > 1, 



(4.31) 



i.e. for |t| ^ 1 the average value ot the position operator behaves classically (see Fig. 7). We 
also point out that the expectation value of the particle velocity {v{t)) = —eT^{^/Xfit) lies 
between -1 and 1, that is it does not exceed the light speed; for |t| ^ 1, |(i^(t))| — > 1, i.e. it 
reaches asymptotically the speed of light. This is yet another evidence of the correctness of 
our approach based on the Salpeter equation. The probability current which can be derived 
with the help of (3.36) is too complex to be reproduced herein. The time development of 
the probablility current obtained from computer simulations is shown in Fig. 8. 



D. Plane wave solutions 



An easy inspection (compare Ref. [22]) shows that the Salpeter Eq. (2.2) with V = 0, i.e. 
in the case of a free particle, possesses plane wave solutions such that 

0(a;,t) = Ce-^(^*-*'-"), (4.32) 

where £ = & + fc c2 and C is a normalization constant. Taking into account (2.8) and 
(3.8) or inserting (4.32) into (3.23) we obtain the following formula on the corresponding 
probability current: 

j = pv, (4.33) 
19 
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FIG. 6. The plot of the probabihty density p{x,t) = \(f){x,t)\^, referring to the wavefunction (4.29) 
of the massless particle in a linear potential. The parameter ji = 1 and A = 1. The classical 
A-shaped dynamics of the maxima of the probability density is easily observed. 

where p = |0|^ = |Cp, and v is the relativistic three-velocity given by 



The relation (4.33) is a natural generalization of the analogous formula on the probabil- 
ity current which is well-known in the nonrelativistic quantum mechanics. We also recall 
that nonrelativistic counterpart of (4.33) is a point of departure for the hydrodynamical 
formulation of quantum mechanics. 

We point out that in the case of the Klein-Gordon equation the probability current for 
plane waves has the same form as (4.33), that is 




(4.34) 



JKG = PkgV. 



(4.35) 



However, the probability density pkg such that 



Pkg = 



(4.36) 



can be either positive or negative depending on the sign of the energy. 
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FIG. 7. The plot of the expectation value of the position operator in the state (4.29), given by Eq. 
(4.30) with n = 1 m"^ and A = 1 versus time (solid line). The dotted line refers to the classical 
trajectory x{t) = — \t\. 




E. Massless particle in three dimensions 



In this section we investigate a free massless quantum particle in three dimensions de- 
scribed by the Salpeter equation 

i» = ^^(x.*), (4.37) 

where we set c = 1. Taking into account the form of the Fourier transformation (4.7) of the 
solution to (4.1) corresponding to the case of a free massless particle in one dimension, one 
can easily guess the following solution to (4.37): 

0(a;,t) = / d^pe^P-'e-^^IPle-'^IPl, (4.38) 

(27r)2 J 

where C is a normalization constant and a > is a parameter. Hence, we get the normalized 
solution to (4.37) which is a plausible three-dimensional generalization of the solution (4.8), 
such that 

, , , (2a) 2 a -I- it , , ^, 

<l>M^^-^y^———^, 4.39 
TT [r'' -I- (a -I- itj^J"' 

where r = \x\. From (4.39) it follows immediately that the probability density is 

Mx.t) = Wx.0r = ^[p3^4^;±5?^^. (4.40) 

The time evolution of the probability density (4.40) is presented in Fig. 9. It appears that, 

in opposition to the case of a free particle in one dimesion (see Fig. 1), the wavefunction 

(4.39) spreads out. Furthermore, using the definition (3.8) one can derive the following 

formula on the probability current: 

/ 32r^t^ + (r^ + f + a^) [3(r^ - + a'y + 12a^t^ - 8rH^] 

^''X 27rW [(r2 -t^ + a2)2 + 4aH^]2 v ■ ) 

3a^ (r + tf + 1 

The plot of the norm of the probability current (4.41) is shown in Fig. 10. As expected the 
norm approaches zero as time passes. 



V. CONCLUSION 

In this work we study the probability current for a quantum spinless relativistic particle 
based on the Salpeter equation as the fundamental equation of the relativistic quantum 
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FIG. 9. The time evolution of the probabihty density (4.40), where [p] = m ^, [r] = m, and a = 1 

m, showing the spreading of the wavefunction (4.39). 
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FIG. 10. The time development of the norm of the probability current (4.41), where [j] = m.~^, 
and a = 1 m. 

mechanics. The introduced probability current and related probability density show a very 
good behavior free of pathologies occuring in the case of the Klein-Gordon equation such as 
the negative probability density or the lack of the massless limit of the probability current. 
Referring to the nonlocality which is often indicated as a grave flaw of the Salpeter equa- 
tion, we have not observed any ill behavior of the discussed general as well as particular 
exact solutions. Quite the opposite, the introduced probability current satisfying Eq. (3.10) 
excludes in view of (3.12) the particle velocities greater than the speed of light. We reahze 
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that the point remains concerning the lack of the manifest covariance of the theory which 
is usually pointed out as the second main disadvantage of the Salpeter equation. However, 
it is our belief that this flaw can be circumvented by using the preferred frame approach 
which has been already successfully applied for the Lorentz covariant localization in quan- 
tum mechanics [38j, relativistic EPR correlations [39j, and Lorentz covariant formulation of 
classical and quantum statistical mechanics [iQl HI] . 
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